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Page 7^ column 1: In equation ( 58 a), the last bracketed expression 

should be corrected as follovs: 


[332 + 6 _ 4 xq (2 + p2)] 


Psige 7^ coliiinn 2: In equation ( 58 b), the last bracketed expression 

should be corrected as follovs: 
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Page 8, column 2: In equation (45a) the last term should be corrected 

as follows: 


+ 2 = 

Page column 2: In equation (45c), the factor 2 preceding the second 

parenthesis should be deleted; that is, the second term within the 
bracket should read 

- (1 - 2xo)(fi + Gi) 

Page 9) column 2; Equation (46a) should be corrected to read as follows: 
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It is pointed out that the foregoing errors have been corrected in a 

subsequent NACA publication (NACA TN 3 O 76 by Nelson, Rainey, and Watkins). 
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BUCKLING OF THIN-WALLED CYLINDER UNDER AXIAL COMPRESSION 

AND INTERNAL PRESSURE ‘ 


Hy Hsu Lo, Harold Crate, and Edward R. Schwartz 


SUMMARY 

An investigation was made of a thin -walled cyclinder under 
axial compresHion and various internal pressures to study the 
e^ect of the internal pressure on the compressive buckling stress of 
the cylinder. A theoretical analysis based on a large-deflection 
theory was also made. The theoretically predicted increase of 
compressive buckling stress due to internal pressure agrees 
fairly well with the experimental results. 

INTRODUCTION 

The buckling of thin-walled cylinders under axial compres- 
sion and lateral pressure has been investigated by Flugge 
(reference 1) who found that the effect of the internal pres- 
sure on the buckling load is negligible. Fliigge’s conclusion 
is in contradiction to the results of a series of tests, made at 
the Langley Aeronautical Laboratory of the NACA, of two 
curved panels under axial compression and various lateral 
pressures. These test results, reported in reference 2, showed 
an appreciable strengthening effect of the lateral pressure 
on the buckling load of the curved panels. The apparent 
discrepancy between these experimental results and the 
prediction by Flugge 's theory made it desirable to investi- 
gate this problem further. Consequently, additional tests 
were made of a cylinder under axial compression and various 
internal pressures for which results are presented herein. 
A theoretical analysis of this problem is also presented 
which differs from that of Flugge in that the present analysis 
is based on large-, rather than small-, deflection theory. 

APPARATUS AND PROCEDURES 

Test specimen. — The specimen used for the tests was a 
cylinder, 32 inches long with a 15-inch inside radius, made 
of 24S-T aluminum alloy sheet of 0.0249-inch average 
thickness. It was closely riveted around two heavy steel 
rings, one at each end. Th(‘ butt joint of the two longi- 
tudinal edges was covered both inside and outside by straps, 
0.032 inch thick and VA inches wide, along the total length 
of the cylinder. (See fig. 1.) 

The two heavy steel rings were made of ){- by 4-inch steel 
bar stock rolled to the diameter of the cylinder. Two %- 
by 2-im4i spreader bat’s wore used to reinforce the ring as 
shown in figurti 1. A ring with a flange, machined fiat, was 
fastened to th(‘. }i- by 4-inch st<Hd ring to providi' an even 
bearing surfa<;e on which a st<M‘l cover plate was fitted. 


Spreader bcfrs 



V Longitudinal strain gages 
X Circumferential strain gages 


Figure 1.— Test specimen and strain-gage positions. 

Three steel blocks were placed on top of the plate. The 
applied compressive load was transmitted from the machine 
head through the three steel blocks to the cover plate. 
The joint between the cylinder and the cover plates was 
sealed. 

Equally spaced along the inside circumference of the cylin- 
der at midlength were 16 strain gages, and directly opposite 
to them on the outside were 16 more gages. These gages 
were placed to measure strains in the longitudinal direction. 
Six more gages, three inside and three oiitsid<\ were placed 
to measure the circumferential st ruins. 

Test procedures. — The specimen was subjected to com- 
pressive load in the 1,200,000-pound universal testing ma- 
chine of the Langley vStructiires Research Laboratory. Com- 
pressed air was used to produce internal pressure, which 
could be ma’-^ ained at any desired constant value. The 
pressure was measured by a manometer. The strains were 
recorded by standard electric strain-gage equipment and 
the end-shortening was measured by dial gages. 

The cylinder was prcloaded and the strain-gage readings 
wove taken. The three steel blocks were so adjusted that 
all longitudinal strain-gage readings around th<'- (dreumfer- 
ence of the cylindtu* were equal. 
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The compressed air was then lot into the rylinder until 
he desired internal pressure was reached. The axial com- 
)ressive load was increased in increments until buckling was | 
observed. At each load increment, all gage readings were | 
ecorded. The load was then decreased until the buckles 
lisappeared and increased a second time to check the reading 
>btained the first time. During all these steps the internal 
pressure was maintained constant. 

The axial load was then reduced and the internal pressure 
was changed to another value. For each value of internal 
pressure the same procedure was repeated. 

EXPERIMENTAL RESULTS 

A typical experimental result is shown in figures 2(a) 
and 2(b) for the case in which the internal pressure was 1% 
psi. In figure 2(a) the compressive load is plotted against 
the strain -gage readings for four different pairs of gages, 
wdthin the range where the load-strain relation is linear. In 
figure 2(b) the strain-gage reading is plotted for all strain 
gages at three compressive loadings close to the buckling 
load. Figure 2(b) indicates that buckling occurs at a com- 
pressive load of 12,700 pounds between strain gages 22 and 
23. (Note the intersection of the curves at two consecutive 
loadings.) A buckle at this location was observed during the 
test. The compressive load at which this phenomenon 
occurs is considered the buckling load. 

vSince the buckling occurs locally and not simultaneously at 
all the gages, the local buckling strain is obtained by dividing 
the buckling load by the slope of the linear portion of the 
load-strain curve (corresponding to the gage at which the 
buckling occurs. The corresponding stress is the buckling 
stress. The buckling stresses for various internal pressures 
were determined in this same way. 

The results are tabulated in table 1 and plotted in figure 3 
in terms of the two nondimensional parameters 

- 

E t 

-Kfy 

where is the buckling stress, p is the internal pressure, 

R is the radius of the cylinder, t is the wall thickness, and 
E is Young’s modulus. Except for the first test correspond- 
ing to ^==0.1028 in which the cylinder had undergone no 
previous buckling, all the tests were carri(Ml out on the cylin- 
der with possible permanent set. 

THEORETICAL RESULTS 

A theoretical analysis for calculating the buckling stress 
of a cylindrical shell under axial compression and internal 
pressure was obtained by a “large-deflection” theory for 



(u) Linear part of load*strain curve for four typical pairs of strain gagi's. 
FiatTRE 2.— Typical experimental result. Internal pressure, j psi. 



(b) Close to buckling loail. 

Kkjlrr 2.— ConcIiKled. 

wliich details are giv(Ui in the a[)i)(Uidix. 'Flu' Iarge-<ietlcctie 
tlieory was first advanced by \ on Karniau and Isif 
(reference^ 3) in tlie study of buckling stn'ss of cylindrie 
shells under axial compression (l)ut without internal pro 
sure). This theory was s»'' .equeiitly impi'oved by Legge 
and Jones (reference 4). In r(‘f(*rence 3 the buckling stre: 
was shown to depend on whether the load was applied by 
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TABLE 1 

BUCKLING STRESSES FOR VARIOl'S INTERNAL PRESSURES 


ExiH*rimentRl 

Theoretical 

P 

7 

A5cr 

P 

7 


0 

0. 19.36 

0 

0 

0.376 

0 

.01715 

.252 

.058 

.02 

.444 

.068 

. 03425 

.277 

.083 

.04 

.480 

.104 

.0514 

.309 

.116 

.06 

.506 

.130 

.0685 

.350 

.156 

.08 

.528 

. 152 

.0856 

.363 

.170 

. 10 

. 547 

. 171 

.1028 

. 407 , 

.213 

. 12 

,565 

. 189 




. 14 1 

. 581 

.205 



Figure 3.— Comparison of theoretical and experimental results of the buckUng stress at 
various internal pressures. 


rigid loading machine or by a dead-weight machine. In 
the present analysis, the loading machine is assumed to be 
rigid. 

The e.xisting procedures for computation of the buckling 
stress by large-deflection theory involve the solution of four 
simultaneous nonlinear equations for each pressure loading. 
The numerical work is quite lengthy. The method used in 
the present study introduces a fifth equation which governs 
the conditions at which the buckling occurs. The fifth 
equation is based on consideration of conservation of energy, 
which is an extension of Tsien’s buckling criterion given in 
reference 5. Although a solution of five simultaneous equa- 
tions is now necessary, the numerical work is actually 



Figure 4— Theoretical and experimental results showing !ht> incrtMiu nt of huckling 
due to internal pressure. 


reduced to a small fraction of that re(|uiretl if the existing 
procedures were used. This reduction in labor is niadt' 
possible through a proper choice of the parametc'rs in the 
equations and the process of the computations. The results 
calculated by the present method are presented in table 1 
and are represented by the solid-line curve in figure 3. The 
curve is cut off at a value of 0,605, corresponding to 
0.169. This constant value of ?cr= 0.605 for ];>0.169 
is the same as that obtained by the classical theory. 

DISCUSSION AND CONCLUSIONS 

From the theoretical and experiimuital results shown in 
figure 3, the internal pressure is seen to have an apprecialile 
strengthening effect on the cylinder. Although the two 
curves obtained from theoretical and experimental results 
do not coincide, both sho>v the same trends as regartls th(‘ 
effect of internal pressure on the buckling stresses. If tlu‘ 
increment of the buckling stress ^cr due to the presence of 
internal pressure (that is, the difference between the buckling 
stress with the pressure 'acr and that without the pressure 
(t^cr)p-o) is plotted against the internal pressure, as shown 
in figure 4, a good agreement is obtained between the theo- 
retical and experimental residts. These data indicate that, 
although further improvement of the theoiy is necessary 
for the determination of the magnitude of the compressive 
buckling stress, the theory gives a fairly good prediction of 
the increase of buckling compressive stress that may be 
expected as a result of internal pressure. The discrepancy 
between the theoretical curve and the experimental curve 
of figure 3 is believed to be caused by such factors as manu- 
facturing imperfections in the specimens, material irregulari- 
ties, and energy absorbed by the loading machine, which 
have not been included in the theory. 

Langley Aeronautical Laboratory, 

Nation/’ Advisory Committee for Aeronautics, 
Langley Field, Va., October 12, 1949. 
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THEORETICAL ANALYSIS OF BUCKLING LOAD OF CYLINDRICAL SHELLS UNDER AXIAL COMPRESSION AND INTERNAL 

PRESSURE BY LARGE-DEFLECTION THEORY 


BACKGROUND OF THEORY 

The use of large-deflection theory for shells under axial 
compression was first advanced by Von Karman and Tsien 
(reference 3) in an attempt to explain the discrepancies 
between the buckling loads predicted by classical theory and 
those obtained from experimental results. (See^ for instance, 
reference 6.) The results of reference 3 indicated that 
cylindrical shells can be maintained in equilibrium in the 
buckled state by a compressive load considerably lower than 
that predicted by classical theory. A plausible explanation 
of this result is that, before the classical buckling load is 
reached during a test, the cylindrical shell ‘‘jumps^^ from an 
equilibrium unbuckled state to an equilibrium buckled state. 
The physical phenomena of the jump were further examined 
in reference 5 by Tsien. 

The treatment of Von Kdrmdn and Tsien in reference 3 
was left incomplete, however, in that the equilibrium posi- 
tions at the buckled state were determined by differentiating 
the total potential energy with respect to some but not all of 
the physical parameters involved. The resulting equations 
gave a relation between the average compressive stress <r and 
the end-shortening e in terms of the remaining parameters. 
A set of curves of a against < were thus obtained for various 
combinations of the remaining parameters. 

Improvement of the theory of Von Kdrmdn and Tsien was 
made by Leggett and Jones (reference 4), who took the 
derivatives of the energy with respect to all the parameters 
and thus obtained a single curve between a and €, represent- 
ing all equilibrium positions of the cylindrical shell in the 
buckleil state. The same result was obtained by Michielsen 
(referen(‘e 7) in a similar process. Such a curve is shown 
by BC of figure 5. 

Theoretically, when the cylinder is compressed, the relation 
between <r and e follows the straight line ODA which repre- 
sents the unbuckled state and will reach the point A if 
everything is perfect; the cylinder then buckles and the 
relationship follows the curve ABC which represents the 
buckled state. Before point A is reached, however, some 
external disturbance may possibly cause the cylinder to jump 
from the unbuckled state represented by the point D to the 
buckled state represented by the point E. The positions 
of D and E on the respective curves depend on the actual 
physical conditions of the jump. 

If the physical condition which governs the jump is known 
or defined, the l)uckling stress corresponding to the point D 
can Ix^ ol)tained directly without going through the labor of 
finding the curve* ABC. This procunlure can greatly reduce 
the amount of numerical work. 



Figure 5.— Relation between the average compressive stress <r and the eiul-shorttMii-u' e. 

In reference 5, Tsien introduced a criterion which gov(‘rns 
the jump DE for the condition of loading obtained in a 
rigid testing machine; namely, that the strain energy remains 
the same before and after the jump and that the jump occurs 
at constant end-shortening. According to this criterion 
the line DE must be vertical and must cut the curve ABC 
in such a way that the two shaded areas ADG and GBE 
are equal. In fact, the area ADG represents the additional 
enei*gy that is needed to assist the cylinder in jumping from a 
condition represented by D to that represented by G and 
the area GBE represents the energy that is given up by the 
cylinder when it arrives at the lower energy level, point E. 
The energy represented by the area ADG is very small, 
and therefore a slight disturbance from the surrounding air 
might assist the cylinder to jump from the unbuckled state 
to the buckled state at a compressive stress well below the 
classi(‘al buckling stress corresponding to point A. 

Since the external disturbance is rec[uir(Ml to assist tin? 
cylinder to jump from the state* (*orn*sponding to D to that 
corT*^sponding to G, a slightly larger external disturbance 
can well cause the cylinder to make the transition from 
the state represented by D' to that represented by B, 
except that in the (*ase in which the cylinder jumps from 
D' to B the cylinder absorbs the energy of the external 
disturbance and does not ro-emit it. ! he buckling stress 
can be then as low as point D^ This fa(‘t was pointed 
out by Tsien in r(*feronce 8. 
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In addition to the two criterions just mentioned, there 
are still othei-s that might be used. In view of the fact 
that the choice of the buckling criterion is a much less 
important factor in the determination of buckling stress 
than are such other factors as, for example, the initial 
imperfections, Tsien^s criterion of reference 5, as represented 
by the line DE, is as reasonable as any other, and the 
choice of this criterion greatly simplifies the numerical work. 

Tsien^s criterion of reference 5 cannot be applied directly 
to the present analysis, however, because with the presence 
of the internal pressure the strain energy is no longer the 
same before and after the jump. In addition, the criterion 
is applied herein in quite a different manner from that of 
reference 5. In reference 5, a series of values of wave num- 
l)er n and aspect ratio were chosen and the criterion was 
applied to each pair of values of n and /8; the pair of values 
of n and /3 which gave the minimum value of buckling load 
was considered t<» correspond to the buckling condition. 
In the present analysis, since the variation of <r with « can 
be plotted only as a single curve, this criterion need be 
applied only once for each internal pressure. The results 
correspond to the minimum-potential-energy condition. 

In the derivation of the present analysis, the basic equa- 
tions in reference 3 are first extended to include the effect 
of internal pressure, Tsien^s criterion governing the jump 
for rigid machine loading (reference 5) is modifit'd, and the 
buckling stress is finally obtained. 


SYMBOLS 


A list of symbols follows. Most of the symbols used in the 
present report are the same as those in reference 3 ; exceptions 
are the use of for Poisson’s ratio, X for wave length, and 
for aspect ratio of the buckled waves. 




/o) f\ , J 2 


V 

t 

y 


u 


length in longitudinal 
length in circumferential 


half wave 
direction 
half wave 
direction 

parameters used in deflection function 
number of waves in longitudinal 
direction within length equal to 
circumference of cylinder 
number of waves in circumference 
internal pressure 
thickness of cylinder wall 
coordinates measured in longitu- 
dinal and circumferential directions, 
n»sp{H‘tively 

component of displacement of a point 
on median surface of shell in 
ar-direction 

component of displacement of a point 
on median surface of shell in radial 
direction 


a 



E 

5 /, B,\ . . . B,') 


E\j E2} • ■ • ^5 




dD 


R 

W, 

W 2 

IT 3 

W, 


measure of average circumferential 
stress per wave length in longitudinal 
direction 

end-shortening of cylinder 
average compressive stress 

aspect ratio of buckled waves 

Poisson’s ratio 
total potential energy 
Young’s modulus 
strain energy 

certain functions of ^ 

certain functions of p and (3 

{D represents the functions Z> 2 , 

. . ♦ Z>,) 

(D represents the functions Di, D 2 , 

• * .• 

radius of cylinder 
elastic extension al energy 
bending energy 

work done by applied compressive load 
work done by internal pressure 


Nondimensional parameters: 


if=n^ 


t 


R 

_ <rR 
"’~E t 


(!) 


w 




<l> 


1 


EtK\t 


(F„Tr,,lF3,and W, are defined in the 
same manner.) 


1 


R 


<r'=(<T— i r,Di 


Subscripts: 

0 

u 

cr 


pertaining to buckling condition 
unbuckled state just prior to buckling 
buckling condition 
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DERIVATION OF BASIC EQUATIONS 

: :vs;;=“' - 




— £’fXaX 


i_ = 4 [( 1 - fj2+^ + 












J 


IV' 




+ 9^')^ (9 + 1 9(1+ 

j/,2 l^i {\+ff^y+\ (l + ^^)]+(l +0^)/./2+(l +^')/^'{ 


n) 


^=[2 (1 (g+n^ I (M+^^) i ^0] 

\ax,x *- 


(3) 


where a, b, fjL, and v have been changed to X^, X^, and m, 
respectively, to agree with the notation of the present 
report and where 


A^l/yn^-Q/.+f.) 


the work done by the internal pressure becomes 

H>=-4;>BX„xy/o+4‘) 

If the total potential energy 

()>=\\ I 11 2 — 11 3 11 p 




and 


<? = ^/l»'(|/l + /2) 

H=ny|/i+/2y 


is differentiated with respect to /o and the derivative is set 
equal to zero, the following expression is obtained: 




(6) 


The work done by the internal pressure is, for a complete 
wave panel, 


Substitute this expression into equations (1), (2), (3), and 
(5) for iri, W., 11^3, and W„, respectively, and the following 
equations are obtained: 


,*(^Bjy+B,jyj2+ Bjyjy+Bj, jy + )- 


\Et\akt 


— rr wp dx dy 


The negative sign is introduced 1 >cause the radial deflection 
w is considered positive inward. 


If the same deflection function given in equation (16) of 
reference 3 is used, that is, 


It 


- =(/»+ {O''' 2 ^ ^ + 4 r) 


-R 


/,/ 2mx . 2ny\ 
•J(eos-,^-rcos-^j 


)+ 

(4) 


W. 




^ + A) I 8 


+ 


ir, 

is: 




n+ 


Ki)«v (+’+■'' 


(9 
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H' 




where 




B I 1 

‘ 64L 8 4 (1+/30“‘^(1+9/3')''^(9 + ^TJ 

R-l p ■ 0* , 

16L2(1+/3 




g__irii I 1 

’ 16 L 2 {\+0T {Q+0yj 

B = - - ^*- 

8{l+pr 

^’=6(r^o [I (1 +^^)] 

1 


B, 


' 6(1 


(l+^‘) 


Equations (7) to (10) may be expressed in terms of the 
nondimensional parameters 7, p, p, t;, f, and W as 

W,^4(7^+2^-5p+y)+r,^^*n,+ v^\-D,)+^W, ( 11 ) 

( 12 ) 

U 3 = 8it"+ 8p<r p+ (X t; (13) 


Wp=Siia p+8p^-^, 


(14) 


whore 


D,-i «■(!+, + ,■) 

D2=5i+^2p +^3P‘ + ^4p'^ + ^ B4P* 

D,==(2B4+~yi+2p) 

Z>4 = ^254 + — ^ + ^ (p + P^) 

D~B^-\- Bfi(p + P*) 

The nondimensional total-potential-energy parameter 0 is 

= -4(P+2m7/5+p^)-(7-1p) ,r=A + 

u*f*£>3-,r’A+m+T,^rz?3 (15) 

The nondimensional strain-energy parameter is 
0 -TFi+TFo 

^4(V^ + 2fiap + y)+r^^D,~ 

(Hi) 


The relation between the end-shortening « and the average 
compressive stress can be determined from equation (23) of 
reference 3 by integrating; thus, 


1 

*=-^ ' 


=-- r 

*/ 0 

^ (e e) l6 (i ^ 

where a/£*, as determined from equation (24) of reference 3, 
together with equation (6) herein is 

a_pB 

E~Et 

Therefore, the relation between 1 and ^ in nondimensional 
form becomes 

- R 


= (T-|-/x/? + - y] f^Z>i 


( 1 “) 


Equations (15), (16), and (17) are for cylinders in the buckl(‘d 
state. For cylinders in the imbiickled state, the correspond- 
ing equations are 

0u= — 4((7u^+2MO’„/j + /r) (IS) 

0^ = 4(^„' + 2/uff„7> + F) n H) 

^M = '^a + /up (20) 

EQUIUBBIUM POSITIONS OF CYLINDERS IN BUCKLED STATE 

The equilibrium positions of cylinders in the buckled state 
can be obtained by differentiating the total potential energy 
of equation (15) with respect to each of the parameters 77, 
f, p, and and by setting the derivatives equal to zero. 
Four simultaneous, nonlinear equations are thus obtained: 

0=— |^(ff— ^ p) riDx-2{ntyO<\- 1 . 0 ( 1 ; s')Z>3 — 

Zy, — 77*L1 ,iJ (2i' ) 

= d = — p^ ’?(Ci). — (?;j')-(Z>j)j — (7)S-)(Zl:,)^ — 

(z>4).-.im)p] r' 


d0 


d(p 

dp 


§|=o=-[(^-^ p) v(D,},-(viy(P,h-(r,n(D,),- 


( 21 ) 
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where 
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(D,),=|^'(l + 2p) 


W,^B,+ 2B,p + ^B,p^+2By 
(Dj),=2 


(D4),=§a+2p) 


(Di),= B,{\+2p) 

(A)p=^*(|+P + '’ 0 =^^' 


(D,)p=B,'+B,'p + B^'p^+B,'p^+^B,'p* 


(P,V=2B4'(l+2p) 
{Dt)p=2B: 
{D,)p=B/+B,'{p + P^ 


and 


- 1^ [?'^T (9+3 J 


r,, 1P__^ + — ^ 
"4 L2 (1 + 3")’^(1 + ®^')' 


3^_ "I 

(^37J 


^ (9+3 J 


1 3* 

®‘'~2(l+30’ 




B.' = ; 




3* 


Let 


r=(<r-ip)rjl>. 


The four simultaneous equations (21) become 

7'=(,,f)’(2 D 2 )— (’ir)D3+2'T^f^5 


7' = (n r )*(2 O 2 ) - (’f f )( 1 ■ 5 ^> 3 ) + D 4 + 

> 


Theoretically these four simultaneous equations <a 
I Tw 1 > « and 3 in terms of 7 for a given pressuiv. It 
^7.^0 b ii ;..«i too o<i.to.oE (iTi. . 
llSsLrloning 7 .od .ho »mP-.ivo s,„,A a |S » 

which represents all equilibrium positions at thi 

state fact, this solution is essentially that obtainec _ 

state in , and Michielsen (reference < ) 

for lyUEdera with «i.l comprrsaion but no inl.ru.l pre.s^ 
PrwticaUv, however, the solution of the four simultaneous 
4 w f221 reauires a long and tedious numerical process. 

k necessary This solution can be obtained b\ tht in 
duction of “one more equation which governs the condition 

at buckling. 

BUCKLING CRITERION 


In relerenoo 5, Taien gave the following criterion wluob 
goverL the condition at buckling: That the .train-energ. 
5 IhTboclded evlinder is the same as the s.tatn energ, oi 
Utekled c'vlinder when the cylinder ts tested ,n tfo 
ritrid^ling machine so that the ettdothortentng Joe. no. 
ciSngc during buckling. This cnlerton is appaienth 

t^tabUshed front consideration, of “"""“7" 

Although other physical criterions can be ustd (for instant 
see reference 8), the criterion of reference .5 was chosen ami 
extended to include the case for which the eternal piessui. 
is present. The choice of this criterion simplifies the numei 

''"^When^nternal pressure is present in the cylinder, work i: 
done by the pressure during buckling. The strain energy 
the buckled state is no longer equal to that in the un uc 

but — ( 

^=^.+MT, <2.5 


where AlT, is the work done by the pressure during bucklinc 


or 


SfXnXa ^ 


ifl Aft 


(22a) 

(22b) 


=Wn-(irnh 

Equation (14) can be rearranged as follows; 


(2 


VTp = 8(p- + Mo-;^-^2 


Therefore, for the unbuckled state, the last term is eliminat 
and 


T = [(, mD 2 )- iv r)(i> 3 )p+( 04 )p+ vW,] (22c) 


7'= [(, f )"(Oa)p-(’J r)(i?3)fl+ (C>4)p+ D5 )p 1 (22d) 


(W’p)n = 8(p- + M<ruP) 

Then, from equations (24), (25), and (26) 


A IT, = Spp(<r 
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'Phe buckling criterion becomes (equations (23) and (27)) 
or, from equations (16) and (19), 

4(7^+ 2m^^+F)+ ^ f 

= 4(7„^+2/i0'„/>-{-;/) + 8(7— ff„)/Ap— i (28) 

Since the end-shortening remains unchanged during buckling^ 
that is, €=Tii, the following relation is obtained from equa- 
tions (17) and (20): 




(29) 


If this relation is substituted in equation (28) and if the 
relation 7'— is used, the buckling criterion 


becomes 


?'=(H)^(£>2-3^W)-(^f)Z)3+04+^^A 


(30) 


The solution of the five equations (22a), (22b), (22c), (22d), 
and (30) gives the buckling stress for a given internal pres- 
sure. The following section presents a very simple method 
for the solution of these five simultaneous equations. 


METHOD OP SOLUTION 


From equations (22b) and (30) and equations (22a) and 
(22b), the following equations are obtained: 


nr 



(31) 


Z>3(nr) 


D, 


(32) 


<^or a preassigned value of /8, assume various values of p and 
ompute ijf and from equations (31) and (32). Substitute 
hese values in equations (22a) and (22c) to obtain (o-')a and 
(r')cy respectively. Plot both ((r')o and (?')c against p. 


The intersection of these two curves determines a pair of 
values <t' and p which are called (x\ and p®. The correspond- 
ing values of {r\^)o and (n^)^ are computed and substituted in 
equation (22d) from which the pressure p can be calcu- 
lated. For each assigned value of there are obtained cor- 
responding values of 7^ and . A curve of against p 
can thus be determined. If the following relations are used. 

the relation between acr and p is obtained as shown in figure 3. 

CUT-OFF BUCKLING STRESS 

When equation (31) is derived from equations (22b) and 
(30), a factor (ijf) = 0 is also obtained. If this relation is 
used instead of equation (31), it can be shown that the 
buckling stress ^cr can never exceed the classical buckling 
stress 0.605 which is independent of pressure. 
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